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Abstract 

In a circle (an S 1 ) with circumference 1 assume m objects distributed pseudo- 
randomly. In the univeral covering manifold R 1 assume the objects replicated ac- 
cordingly, and take an interval L > 1. In this interval, make the normalized his- 
togram of the pair separations which are not an integer. The theoretical (expected) 
such histogram is obtained in this report, as well as its difference to a similar his- 
togram for non-replicated objects. The whole study is of interest for the cosmic 
crystallography. 

1 Introduction 

Cosmic crystallography (CC) is a method to unveil the topology of the universe, and 
initially looked for spikes in a pair separation histogram (PSH) Since spikes are 
absent in hyperbolic spaces, it appeared that the method was useless in such spaces. 
However, it was soon shown that not only a Clifford translation (responsible for a spike) 
press its fingerprint on a PSH, but also the other isometries of the space 0. 

When spikes are absent, the PSH of a ball containing repeated images - the (p m (l) 
- is very similar to that of a ball with same radius and same geometry, but without 
duplication of images - the 4> S {1). A suggestion was then made, of studying the difference 
of the multiply and the simply connected histograms, <fi m (l) — (fi s (l) 0. 

To improve the method, expected functions 4> s exp {l) were derived to replace the his- 
tograms 4> s (l) obtained from computer simulations, for all three geometries with constant 
curvature Q. Graphs of (p m {l) —(f> s exp (l) were obtained, clearly evincing the topology of an 
euclidian, an elliptic, and a hyperbolic three-space ||. The contribution of each individual 
isometry g to a PSH was examined, and normalized histograms <p 9 (l) (denned in ref.0) 
were obtained from computer simulations ||; these simulations also gave histograms of 
4> U {1) — 4> s exp (l), a previously unsuspected quantity ||. 
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Recently the exact (noiseless) functions (p 9 exp {l) were given for the euclidian isometries 
0. In the present report we finally have a first acquaintance with functions 4>e Xp (l), the 
exact (noiseless) counterparts of the 'uncorrelated' normalized histograms <p u (l) defined 
in ||. We examine a one-dimensional system: a universe with topology S 1 , a circle with 
circumference 1; we assume the horizon at a distance L/2 on each side of an observer, so 
the visible universe has total length L; clearly if L > 1 then there are repeated images 
in this visible universe. In section 2 we give a detailed description of how to obtain the 
expected uncorrelated signature (fL exp (l) when 1 < L < 2. In section 3 we exhibit the 
generalization for arbitrary horizon L/2. In the Conclusion we make a few comments, 
and in four Appendices we derive a few somehow lengthy mathematical results stated in 
the report. 



2 When 1< L < 2 

In a computer simulation, we usually execute the following set of prescriptions to obtain 
the uncorrelated signature <Pl{1): 

m - p 

P / \\\ P 



Ox 1 L 

Figure 1 The distribution of objects in the interval (1,L) is an exact copy of the distribution 
in (0, x); here p = 3 and m = 8. 



1. in an interval (0, 1) randomly distribute m objects; see figure 1; 

2. in the side interval (1, L) make an exact replica of the p objects laying in (0, x); 

3. measure the (m + p) (m + p — l)/2 separations I between the total m + p objects, 
and discard the p correlated separations (those which have I = 1 exactly); 

4. make a normalized histogram of the 

V mp = ~(m + p)(m+p- 1) -p (1<L<2) (1) 

uncorrelated separations; 

5. make a large number of new normalized histograms, by repeating the steps 1 to 4 
with same m (although p usually varies); 

6. take the mean of these histograms, < 4>mL^) an d construct the quantity 

< cp- mL (l) >= ( n - 1 - £ u a ) [< Cl(0 > -Ml)} , (2) 
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where 



,(() 



id-!- 



< / < L, 



(3) 



and where the factor n — 1 — Y^ u g — ( m ~ 1)-^ ~~ x 0- ~ X )/L is explained in the 
appendix 1; 

7. the (computer simulated) -uncorrelated signature < ^(0 > is the quantity < 
Vmi(0 > when m — > oo; in practice m > 50 usually suffices. See figure 2. 





Figure 2 Computer simulated functions < ^PmL^) > f° r { m = 2,L = 1.7} and {m = 30, L = 
1.3}. 

We now develop an analytical method to obtain the uncorrelated signature v^lIO- We 
are dropping the subscript exp in all expected (theoretic, analytic, mean) probability 
distributions. Initially define the lengths x and y (see figure 1) 



x = L-l, y = l- x (1<L<2), 



(4) 



and assume that m objects are randomly distributed in (0, 1); the probability that p 
objects be in the interval (0, x) and m — p objects be in the interval (x, 1) clearly is 



V 



mpx ^ m 



u m x y , u 



ml 



p\{m —p)V 



irrespective of the values of m and x we have 

771 

^ t "Pmpx ~ 
p=0 



(5) 



(6) 



We denote as <^ pL {l)dl the probability of finding in (0, L) an uncorrelated pair with 
separation between I and I + dl, when there are m objects in (0, 1) and p objects in (0, x); 
clearly it satisfies 

rL F mpL {l)dl = l. (7) 

Recall that a pair (P, Q) is said (^-correlated when the isometry g brings one of the 
members to the other; the pair is uncorrelated when no such g exists. To investigate 
0" L (l) when 1 < L < 2 we first call A the interval (0,x), call B = (x, 1), and call 
C = (1,L), and note that there are 
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• waa = p{p — l)/2 pairs with both members in A; 

• wab = p{fn — p) pairs with a member in A and the other in B; 

• wac — p{p — 1) uncorrelated pairs, with a member in A and the other in C; 

• Wbb — { m — p){itl — p — l)/2 pairs with both members in B; 

• wbc{= wab) pairs with a member in B and the other in C; 

• wcc(= waa) pairs with both members in C . 

In total, there are V mp (eq. (]!])) pair separations to be considered. 

A short reflection gives that the density 4>m P L^) can be decomposed as 



r mpL {i) 



V 



mp 



waa4>aa{1) + wab4>ab{1) +wac4>ac{1) + 
w B b4>bb{1) + Wbc4>bc{1) +wcc<ficc(l) 



where each 4>xy(1) is the probability density of finding an uncorrelated pair of objects 
separated by /, one in X and the other in Y; clearly all obey 



<j> XY {l)dl = 1. 



(9) 



There are two basic types of according as X = Y or X ^ Y. When X — Y, 

suppose a segment of length fi, and randomly select two points of it; the probability that 
their separation lie between / and I + dl is cj) s ^{l)dl with (see figure 3) 



< / < ji. 



(10) 



2/n 




n I 

Figure 3 Pair separation density function for an interval \x. The underlying area is 1. 

When I / 7, consider two intervals with lengths a and (3, with separation 5 (see 
figure 4); randomly select one point in each a and (3; the probability that the separation 
between these points lie between / and I + dl is 4>s{a/3)(l)dl, with the density <fis(ap){l) as 
depicted in figure 5. 
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Figure 4 Intervals with lengths a and j3, with separation 5; assume a < (3. 

Figure 5 The probability density <ps(a^){f) for a < /3 (see figure 4); three particular cases are 

also displayed; all underlying areas are = 1. 



The functions 0xy(O appearing in eq.fll) are as displayed in the figure 6, for the case 
with x < y; for x > y a similar set has to be constructed, see figure 7. 




When x < 0.5 the density 0^1,(0 > ec l-®> i s a sequence of four straight segments with 
endpoints at I = 0,x,y, 1, and L (in this order), and values 

C P x(0) = ^-[w BB ^ + 2w AA h (11) 
C p lW = ^-[wbb 2 ^~ X) +2w AB ~] (x<0.5), 
r mpL (y) = ^-[2w AB k (z<0.5), 

L^mp y 

C P x(l) = J-[v>Ach <p u mpL (L) = 0. 

L^mp X 
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When x > 0.5 the sequence of endpoints changes to I = 0,y,x,l, and L, and the 
values of 0^pl(O at / = ?/ and I = x become 



r mpL (v) 



mpL 



1 [2w AA 2{x ^ V) +2w AB ±] (x>0.5) 



x 



V L ~ M x 2 

1 1 x — y 

— — [2w AB - + w AC — — 



(12) 



(x > 0.5). 



Two examples of functions f° r ^ < L < 2 are shown in figure 8. 




x y 




y x 



Figure 8 The probability density ^pl(0 f° r m = 3,p = 2, and two values of L: 1.4 and 1.6 . 
Both underlying areas are 1. 

Having the m + 1 functions 4>m P L^)y P = 0> • we introduce the probability density 

m 

Cl(0 = E^C p l(0, (13) 

whose interpretation is obvious: (j}^ L {t)dl is the probability that two uncorrelated objects 
randomly selected in L have separation between I and l+dl, when m objects were randomly 
distributed in the interval (0, 1). Examples of </Cl(0 are given in figure 9. 



m= 3 
x— .6 




z 



Figure 9 Probability densities <?^l(0 for 1 < L < 2. The graph of 4> s L (l) is given in dotted line, 
for comparison. 

Cosmic crystallography is mostly interested in systems with m » 1. In this limit the 
function closely resembles the simple triangular function 4> S L {1) (eq.@, fig. 3), so 

one is led to define the difference 



lim mL 

m— >oo 



Cl(O-01(O 



(14) 



the asymptotic uncorrelated signature of L. 
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We soon find that the function (p^(l) has a number of symmetries: 

<fl(0) = <fl(L/2) = <fl(L) = 0, <ft(x) = -<fl(l). (15) 

In other words, every <^l(0 with 1 < L < 2 is composed of three line segments, with the 
first segment parallel to the third (see figure 10). As expected, the entire graph of ^(0 
is uniquely fixed by the number f(L), the value of ^(Z) at I = x; in the appendix 2 we 
show that 

f{L) = ^ (1<L<2). (16) 
A plot of f(L) valid for arbitrary L > 1 is given in figure 11. 




Figure 10 Geometro-topological signature (f^(l) for L =1.1 , 1.5, and 1.9 . 

Figure 11 The function f(L), the absolute maximum of ^(Z) (which occurs at I = x); three 

particular values of L are marked, those used in figure 10. 

3 When L > 2 

The generalization of the previous results for arbitrary values of L is straightforward but 
lengthy, so we only state the final results in this section. See the appendix 3 for details. 
The graph of the uncorrelated signature ( |14D with 

L = X + x, XeZ + , 0<x<l (17) 

has the aspect of a slanted saw; see figure 12, drawn for A = 5 and x = 0.2. 
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There are A maxima, which occur in the positions I = x, 1 + x, (A — 1) + x , and 
there are A minima, which lay in the positions I = 1, 2, A. A straight line connects the 
maxima, another one connects the minima, both have angular coefficient —8xy/L 3 . The 
A + 1 segments with positive angular coefficient are parallel, as well as the A segments 
with negative slope. As expected, the value of L is the sufficient datum to draw ^(1), 
since 

<pl(x) = -<pl(\) = (18) 
as shown in the appendix 4. The graph of f(L) = 8Xxy/L 3 is given in Figure 11. 



4 Conclusion 

In our first contact with the cosmic crystallography it appeared plausible that the nor- 
malized expected functions 4>e Xp (l) and 4> s exp (l) were the same, since both are concerned 
with separations between objects isometrically unrelated ||. However, in our computer 
simulations a persistent non-nullity of the difference < <fi u {l) > —4> s exp {l) made imperative 
a more close exam. It soon became evident that a difference indeed existed, and that it 
diminished as the number n of objects present in the sample increased. 
Further investigation suggested to define the uncorrelated signature 

V u exp (l) = (n - 1 - £ v B ) [€ xp (l) ~ € xp il)] , (19) 

where v g = N g /n , with N g =number of g-pairs in the observed universe; for the cosmic 
crystallography we usually have n » 1 + J2 v g- 

Earlier attempts to find 4>e Xp {l) for three-dimensional balls failed, and also for 2D 
balls; we then focussed our attention on a ID ball, this report. When we compare the 
final theoretical result ( |bj|) with the mean of an increasing number of histograms obtained 
from computer simulations, we note a rapid agreement of the two approaches in the region 
of large separations I > L/2, while in the region where I < L/2 a quite larger number of 
simulated catalogs is demanded. This can be seen in Figure 2, where we observe that the 
statistical fluctuations for I large are sensibly less pronounced than those for small /. 

When L < 1, then there is no replication of objects; in this case 0^(/) = 01(0 and 
clearly = 0. When L > 1 is an integer, then objects are replicated; nevertheless still 
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0l(O = 01(0 an d = 0- This can be seen in the figure 11, where we note that f(L) 

vanishes for L = integer > 0. 

Appendix 1 

We evaluate the quantity n — 1 — J2 v g for a universe S 1 with circumference 1 and 
observed universe with total amplitude L = X+x, being A a positive integer and < x < 1. 

Assuming m objects along the circle S 1 with radius 1/(2tt), then the expected number 
of objects in L is n = mL. The sum v g — v -\ + ^-a+i + ••• + v~\ + V\ + ... + ^a-i + ^a 
indeed simplifies to 2{v\ + f 2 + ... + fx), since z/_j = i^. 

Now remember that for z a positive integer is the expected number of pairs of 
objects in the observed universe whose separation is % < A [0; its value is 

nui = m(L — i) . (20) 
As a consequence v g — ^(L — y) / L, and finally 

n-l-Y,v g = (m-l)L-^. (21) 

Appendix 2 

We show that f^x) = 8xy/L 3 when 1 < L < 2: from (11) or flT2|) we have at I = 1 

«WU = (22) 

so we have from 



fc(l) = ~E^(p-l), (23) 

whose value is sought, correct to order m _1 when m » 1. In this limit we have 

E P mp ,(p/m) fc = a; fc + * ( *~ ^ ys*" 1 + 0(m~ 2 ), (24) 

and consequently 

E V mpx F{p/m) = F{x) + ~^F{x) + 0(m- 2 ). (25) 
For m » 1 in eq.(|l]) we find that 

pep- 1 ) 2 ^ 2 i mt+m-i) 0(m _ 2) e ._ , (26) 
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so from ( P5| ) we obtain 
p(p- !) 



p=0 



2x^ 



2x(2x + l)(x - 1) xy d 2 



mpx 



V 



nip 



1 + x 



mil + x) 4 



2m dx 2 



2x 2 

T+x) 



0(m 



2x 2 8x 2 y _ 2 



Since 01,(1) = 2x/L 2 , we finally have from (|l4] 

8xy 



L 3 



>J), and (12 
1 < L < 2) 



(27) 



(28) 



Appendix 3 

We generalize for arbitrary L > 1 the results obtained for 1 < L < 2, in particular the 
equations ([]]) and (0). We first decompose the total interval (0, L) into 2A + 1 subintervals 
according to figure 13, drawn for A = 5. 

A 1 B l A 2 B 2 A 3 B^A^B^Ax \ A Ul 



0x1 2 3 ^-1 XL 

Figure 13 The one-dimensional observed universe with length L = A + x, partitioned into A + 1 
intervals Ai with length x and A intervals Bi measuring y = 1 — x. 

For m objects randomly distributed in the universe (0, 1) we expect p = mx objects 
in each interval Ai and m — p = my objects in each Bi. The number of objects in 
the observed universe (0, L) being mX + p, the total number of pairs of objects in it is 
(mX + p)(mX + p — l)/2; if we deduct the pX(X + l)/2 correlated pairs with members in 
the A's, and the (m — p)(X — l)A/2 correlated pairs with members in the B's, then we 
obtain the expected number of uncorrelated separations (cf eq.(P): 



V mp = ^(mX+ p){mX +p-l)- ^A(A + l)p - ^A(A - l)(m - p). ( 20) 



We next note in (0, L) the existence of 

• w A t Ai — p(p — l)/2 pairs with both members in A{, 

• WB t Bi = (tu — p)(m — p — l)/2 pairs with both members in 

• WA t Aj >t = 2wAiA t uncorrelated pairs, with a member in Ai and the other in A/>j; 

• WBiB^i = ^WBiBi uncorrelated pairs, with a member in Bi and the other in -Bj>j; 

• WA i B J>i = p{™ — p) pairs, with a member in Ai and the other in Bj>i\ 

• WBiA^i = WAiB^i pairs, with a member in Bi and the other in Aj>j. 
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There is a total of (2A + 1)(A + 1) such numbers wxy, and their sum clearly is the V mp 
given in (29). 

With the probability densities 4>xy{1) defined as before, the normalized probability 
density 4>m P L^) * s written similarly to eq. (8), 



mpL 



(0 



v 



w xy4>xy{1)- 



(30) 



mp X V 



As a matter of fact, there are only three essentially different wxy, which we dub waa, wab, 
and wbb, as in sec. 2. Also, there are indeed only 3A + 1 different functions 0xy(O; e& ch 
appearing with variable multiplicity rrixY- These functions, together with the correspond- 
ing mxY and weights wxy, are displayed in figure 14, drawn for L = 5.2. 



W AA 



2/x 
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1+1 21 2(1-1) 6 4 2 

A A A A A 



x y 1 



2 3 
6 



1-1 



y vl 21 2(1-1) 6 4 2 




Figure 14 The 3A + 1 different functions 0xy(O when x < 0.5. On top of each function the 
corresponding multiplicity mxY is written. On the left side the corresponding weight wxy is 
also given. The value L = 5.2 was taken for definiteness. 

When x > 0.5 the set of functions 0xy(O nas a different aspect; see figure 15, drawn 
for L = 5.8. 
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Figure 15 The 3A + 1 different functions 0^y(O when x > 0.5. The multiplicities mjy and 
weights wxy are indicated as in figure 14. The value L = 5.8 was taken for definiteness. 



It is now clear that the functions $L>l(0 are a sequence of 3 A + 1 segments, each 
segment having endpoints either at an integer or separated x from an integer; as a con- 
sequence, also the functions 4>m L (l) (eq.fllBD) have that behavior, as well as the functions 
if u mL (l) (eq.©). See figure 16. 
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x yi 



Figure 16 The function V^l© for m = 2 and L = 5.2. A straight line connects the points with 
abscissa I = i + x (i = 0, A); another, parallel, connects those with I = i + y (i = 0, A — 1); 
also the points with I =integer are aligned. 

Appendix 4 

We generalize eq.(^) for arbitrary L > 1. For m » 1 and T> mv as in eq.(^) we have 
P(P-1) 2£ 2 2Ag(2g + A)(£-l) ^ 

^7 = KW + m ( e + A)4 + ° ( ™ >■ (31) 

while eq. (|27|) now reads 

„, 2x 2 8Xx 2 y ^, , , 

Finally ( [28| ) becomes 

V2(A) = L > 1. (33) 

The graph of /(£) = 8Xxy/L 3 is given in figure 11. 
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